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Abstract 

A path integral approach to quantum physics has been developed. 
Fractional path integrals over the paths of the Levy flights are de- 
fined. It is shown that if the fractality of the Brownian trajectories 
leads to standard quantum and statistical mechanics, then the frac- 
tality of the Levy paths leads to fractional quantum mechanics and 
fractional statistical mechanics. The fractional quantum and statis- 
tical mechanics have been developed via our fractional path integral 
approach. A fractional generalization of the Schrodinger equation has 
been found. A relationship between the energy and the momentum of 
the nonrelativistic quantum-mechanical particle has been established. 
The equation for the fractional plane wave function has been obtained. 
We have derived a free particle quantum-mechanical kernel using Fox's 
H function. A fractional generalization of the Heisenberg uncertainty 
relation has been established. Fractional statistical mechanics has 
been developed via the path integral approach. A fractional general- 
ization of the motion equation for the density matrix has been found. 
The density matrix of a free particle has been expressed in terms 
of the Fox's H function. We also discuss the relationships between 
fractional and the well-known Feynman path integral approaches to 
quantum and statistical mechanics. 
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1 Introduction 



The term "fractal" was introduced into scientists' lexicon by Mandelbrot [T]. 
Historically, the first example of the fractional physical object was the Brow- 
nian motion, whose trajectories (paths) are nondifferentiable, self-similar 
curves that have a fractal dimension that is different from its topological 
dimension [1], [2]. In quantum physics the first successful attempt to apply 
the fractality concept was the Feynman path integral approach to quantum 
mechanics. Feynman and Hibbs [3| reformulated the nonrelativistic quantum 
mechanics as a path integral over the Brownian paths. Thus the Feynman- 
Hibbs fractional background leads to standard (nonfractional) quantum me- 
chanics. 

We develop an extension of a fractality concept in quantum physics. That 
is, we construct a fractional path integral and formulate the fractional quan- 
tum mechanics ^ as a path integral over the paths of the Levy flights. 

The Levy stochastic process is a natural generalization of the Brownian 
motion or the Wiener stochastic process [5], [H]. The foundation for this 
generalization is the theory of stable probability distributions developed by 
Levy The most fundamental property of the Levy distributions is the 
stability in respect to addition, in accordance with the generalized central 
limit theorem. Thus, from the probability theory point of view, the stable 
probability law is a generalization of the well-known Gaussian law. The Levy 
processes are characterized by the Levy index a, < a < 2. At a = 2 we 
have the Gaussian process or the process of the Brownian motion. Let us 
note that the Levy process is widely used to model a variety of processes, 
such as turbulence [S], chaotic dynamics [H], plasma physics [lU], financial 
dynamics [11], biology and physiology [12]. 

As is well known, in the Gaussian case the path integral approach to 
quantum mechanics allows one to reproduce the Schrodinger equation for 
the wave function. In the general case we derive the fractional generalization 
of the Schrodinger equation [see Eq. fl28l) ]. The fractional generalization of 
the Schrodinger equation includes the derivative of order a instead of the 
second {a = 2) order derivative in the standard Schrodinger equation. This 
is one of the reasons for the term "fractional quantum mechanics" (FQM). 

The paper is organized as follows. In Sec. II we describe two fractals: (i) 
a trajectory of the Brownian motion, and (ii) a trajectory of the Levy flight. 
In Sec. Ill we define the fractional path integrals in the coordinate and phase 
space representations. We develop the FQM via a path integral. It is shown 
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in what way the FQM includes the standard one. We derive a free particle 
fractional quantum-mechanical propagator using the Fox's H function. The 
fractional dispersion relation between the energy and the momentum of the 
nonrelativistic fractional quantum mechanical particle is established. 

In Sec. IV the fractional generalization of the Schrodinger equation in 
terms of the quantum Riesz fractional derivative is obtained. The fractional 
Hamilton operator is defined and its hermiticy is proven. 

As a physical application of the developed fractional quantum mechanics, 
a free particle quantum dynamics is studied in Sec.V. We introduce the Levy 
wave packet, which is a fractional generalization of the well-known Gaus- 
sian wave packet. Quantum- mechanical probability densities in space and 
momentum representations are derived. The fractional uncertainty relation 
is established. This uncertainty relation can be considered as a fractional 
generalization of the Heisenberg uncertainty relation. 

In Sec. VI we develop the fractional statistical mechanics (FSM) by means 
of the fractional path integral approach. The main point is go from imaginary 
time (in the framework of the quantum-mechanical fractional path integral 
consideration) to "inverse temperature" it — >■ hl3, where (3 = l/ksT, is 
Boltzmann's constant, h is Planck's constant and T is the temperature. We 
have found an equation for the partition function of the fractional statisti- 
cal system. The fractional density matrix for a free particle is expressed in 
analytical form in terms of the Fox's H function. We have derived the new 
fractional differential equation [see Eq. fl63l) ] for the fractional density ma- 
trix. In the conclusion, we discuss the relationships between the fractional 
approach and the well-known Feynman path integral approach to quantum 
and statistical mechanics. 

2 Fractals 

The relation between fractals and quantum (or statistical) mechanics is eas- 
ily observed in the framework of the Feynman path integral formulation [3j. 
The background of the Feynman approach to quantum mechanics is a path 
integral over the Brownian paths. The Brownian motion was historically the 
first example of the fractal in physics. Brownian paths, are nondifferentiable, 
self-similar curves whose fractal dimension is different from its topological di- 
mension. Let us explain the fractal dimension with two examples of fractals: 
(i) the Brownian path and (ii) the trajectory of the Levy flight. 
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(i) A mathematical model of the Brownian motion is the Wiener stochas- 
tic process x(t) [5]. The probability density pw{xt\xotQ) that a stochastic 
process x{t), will be found at x at time t under the condition that it starting 
at t = to from a;(to) = xq, satisfies the diffusion equation 



dpw{xt\xoto) ^ |v^PH/(a^^ko^o), Pw{xt\xoto) = S{x - xq), = ^ 

the solution of which has the form 

Pw{xt\xoto) = Pw{x -xo;t- to) = = exp{-— — — — }, (1) 

y'2iTa{t - to) 2a{t-to) 

where a is the diffusion coefficient, and t > to. 
Equation ([1]) implies that 

{x-XoY (xa{t-to). (2) 

This scaling relation between a length increment of the Wiener process 
Ax = X — Xo and a time increment At = t — to allows one to find the fractal 
dimension of the Brownian path. Let us consider the length of the diffusion 
path between two given space-time points. We divide the given time interval 
T into slices, such as T = A^At. Then the space length of the diffusion 
path is 

L = NAx = ^Ax = aT{AxY\ (3) 

where the scaling relation [Eq.([2])] was taken into account. The fractal 
dimension tells us about the length of the path when space resolution goes 
to zero. Ax — > 0. The fractional dimension dfractai may be introduced by [1], 



L oc {AxY~'^'''^'''^\ 

where Ax — ^ 0. Letting Ax — in the Eq.([3]), and comparing with the 
definition of the fractal dimension dfractai, yields 

■^{Brownian) _ ^ / a\ 

^fractal ~ ^- \^) 

Thus the fractal dimension of the Brownian path is 2. 
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(ii) Another example of a fractal is the random process of the Levy "flight" 
(or the Levy motion). As discussed in Sec. I, the Levy motion is a so-called 
a-stable random process, and may be considered as a generalization of the 
Brownian motion. The a-stable distribution is formed under the influence of 
the sum of a large number of independent random factors. The probability 
density PL{xt\xQtQ) of the Levy a-stable distribution has the form 



PL{xt\xoto) 



oo 

1 / dfce^'=(^-^o)exp{-a„|fcr(t-to)}, (5) 

2lT J 



where a is the Levy index < a < 2, and a a is the generalized diffusion 
coefficient with the "physical" dimension [cTq,] =cm"sec~^. The a-stable 
distribution with < a < 2 possesses finite moments of order fi, /j, < a, 
but infinite moments for higher order. Note that the Gaussian probability 
distribution is also a stable one (a = 2) and it possesses moments of all 
orders. 

We will further study a fractional quantum and statistical mechanics, 
and it seems reasonable to suggest that there exist moments of first order or 
physical averages (for example, an average momentum or space coordinate of 
quantum particle; see Secs.V and VI). The requirement for the first moment's 
existence gives the restriction, 1 < a < 2. 

The a-stable Levy distribution defined by Eq.([5]) satisfies the fractional 
diffusion equation 

dpL{xt\xoto) ™ / ,| , X ™ - /a\ 
g-^ = pL{xt\xoto), V =-^, (6) 

PL{xt\xoto) = 6{x - Xq), 

where V" is the fractional Riesz derivative defined through its Fourier 
transform [13], [H] 

oo 

Vy{x,t) = -— [ dke'''''\k\''p{k,t). (7) 
27r J 

— oo 

Here p{x, t) and p{k, t) are related to each other by the Fourier transforms 
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p{x,t) = — I dke'-^''p{k,t), p[k,t) = I dxe ^^''p{x,t). 
2n 



Equation ([5]) implies that 

(x-xo) oc (cT^(t-to))^^", 1 < a < 2. (8) 

This scahng relation between a length increment of the Levy process 
Ax = X — Xo and a time increment At = t — to, allows one to find the 
fractal dimension of a trajectory of a Levy path. Let us consider the length 
of the Levy path between two given space-time points. Dividing the given 
time interval T into slices, such as T = NAt, and taking into account the 
scaling relation [Eq.([S])], we have 

L = NAx = —Ax = DTiAxV-". 

At ^ ' 

Letting Ax 0, and comparing with the definition of the fractal dimen- 
sion dfractai P, 0, yields 

diS = «, 1<«<2. (9) 
Thus the fractal dimension of the considered Levy path is a. 



3 Fractional path integral 

If a particle at an initial time ta starts from the point Xa and goes to a 
final point xy, at time t^, we will say simply that the particle goes from a to 
h and its trajectory (path) x(t) will have the property that xita) = Xa and 
x(tfe) = Xh- In quantum mechanics, then, we will have a quantum-mechanical 
amplitude, often called a kernel, which we may write KF{xhtb\xata), which 
we use to get from the point a to the point h. This will be the sum over all of 
the trajectories that go between that end points, and of a contribution from 
each. If we have a quantum particle moving in the potential V{x) then the 
quantum- mechanical amplitude Kp{xijti,\xata) may be written as [3] 
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x{tt)=Xf, 



KpiXbhlXata) = J VFeynmanXir) ■ exp{-^ j dTV{x{T))}, (10) 

x{ta)=Xa ta 

where V{x{t)) is the potential energy as a functional of a particle path 
x{t), and the Feynman path integral measure is defined as 



x{ti,)=x, 



V FeynmanXir).... = \im^ / dXi...dXN-l ( ^ ) X (11) 



N 



nexp|^(x, -x,_i)2| 



X 

here m is the mass of the quantum mechanical particle, h is the Planck's 
constant, xq = Xa, xn = Xb and e = {tb — ta)/N. The Feynman path integral 
measure is generated by the process of the Brownian motion. Indeed, Eq. (fTTil 
implies 

/ \ 

(x,-a;,_i)oc(^-j (At)V2. 

This is the typical relation between the space displacement and the time 
scale for the Brownian path. This scaling relation between a length increment 
{xj — Xj-i) and a time increment At implies that the fractal dimension of 
the Feynman's path is d^^^^^j™'""'' = 2. As is well known, the definition 
given by Eq. lfTT]) leads to standard quantum mechanics. We conclude that 
the Feynman-Hibbs fractional background leads to standard (nonfractional) 
quantum mechanics |3j. 

We propose the fractional quantum mechanics based on the new fractional 
path integral 

x{tb)=Xt t(, 

KL{xbtb\xata) = j Vx{t) ■ exp{- ^ j drV {x{t))} , (12) 

x{ta)=Xa ta 
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where V{x{t)) is the potential energy as a functional of the Levy particle 
path, and the fractional path integral measure is defined as 



x{ti,)=xi, 



J Vx{t)... = (13) 



^ (^a ) — 

-N/a N r 1 / ft \ V" 



where is the generalized "fractional quantum diffusion coefficient" , the 
physical dimension of which is [Dq,] =erg^~°-cm°-sec~", h denotes Planck's 
constant, Xq = Xa, xn = Xb, e = (4 — ta)/N, and the Levy distribution 
function La is expressed in terms of Fox's H function [15] - [17] 




x\ } = (14) 



1,1/a), (1,1/2) 
(1,1), (1,1/2) 

Here a is the Levy index and, as it was mentioned in Sec. II, we consider 
the case when 1 < a < 2. 

The functional measure defined by Eq. f|T3|) is generated by the Levy flights 
stochastic process. We find from Eq. f|T3|) that the scaling relation between a 
length increment {xj — Xj_i) and a time increment At has a fractional form 



\xj - Xj^il oc [h'^-'DaY^" (At)^/". 

This scaling relation implies that the fractal dimension of the Levy path 
is d[^^3^ = a. So, in the general case 1 < a < 2 Levy fractional background 
leads to fractional quantum mechanics. Equations f|T2l) - f|T4|) define the new 
fractional quantum mechanics via the fractional path integral. 

As a physical application of the developed fractional path integral ap- 
proach let us calculate the free particle kernel K^^\xbtb\xata), and compare 
it with the Feynman free particle kernel Kp\xbtb\xata)- For the free particle 
V{x) = 0, and Eqs.([l2]) and ([13]) yield 
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Kf\x^h\Xata) = j Vx{t)-1= (15) 



^ ) " [ ^ \iDa{tb - ta 

It is known that at a = 2 the Levy distribution is transformed to a Gaus- 
sian, and the Levy flights process is transformed to the process of Brownian 
motion. Equation f|T5l) . in accordance with the definition given by Eq. (ll4p 
and the properties of the Fox's function at a = 2 (see Refs.|T6], [17]) is 
transformed to a Feynman free particle kernel (see Eq.(3-3)) of Ref.[3]) 



K',\xM^ata)=^ J -^^^[-Wh^]- ^''^ 

Thus the general fractional [Eq. ffTSj) ] includes, as a particular, Gaussian 
case at a = 2, the Feynman propagator. 

In terms of a Fourier integral (momentum representation), the fractional 
kernel Kf\x\)ti,\xata) is written as 



00 



Kl'{xbtb\xata) = / dp-expU i y (17) 



while Eq. (]T6l) in the momentum representation has the form 



00 



(0)/ +1 1 fj j.piXb-Xa) .piU-ta) 



K, {xbhWQ = ^ dp.exp^ ^ \ ■ (18) 



We see from Eq.f ll7p that the energy Ep of the fractional quantum me- 
chanical particle with momentum p is given by 

Ep = Z^^br. (19) 
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This is a dispersion relation for the nonrelativistic fractional quantum- 
mechanical particle. The comparison of the Eqs. fll7l) and fllSp allows to con- 
clude that at a = 2 we should put D2 = l/2m. Then Eq. ffTOj) is transformed 
to the standard nonrelativistic equation Ep = /2m. 

Using Eq. (1171) we can define the fractional functional measure in the phase 
space representation by 

Dx(r) / Dp(r)... = (20) 



3^(^a ) — 



00 



-00 —00 



, -PNiXb- Xn-i) .Da\PN\"£ 

X exp < I t 

n n 



here e = (4 — ta)/N. Then the kernel KL{xbtb\xata) defined by Eq. (fT2 
can be written as 

00 00 
Kiixhtblxata) = lim / dxi...dxN-i ,^ / dpi...dpNX 

N^oo J [ZTin)" J 



In the continuum limit N ^ 00, e ^ 0, we have 

x{ti,)=x,, r % 

KL{xbtb\xata) = j Dx{t) j Dp(r) exp < ^ j dT[p{T)x{T) - Ha{piT),x{T)] 

x{ta)=Xa V ia 

(21) 

x{t,,)=Xb 

where the phase space path integral / Dx(r) jDp(r)... is given by 

X (ia )— 

Eq. (l20l) . X denotes the time derivative, is the fractional Hamiltonian 
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H^{p,x) = DM'' + V{x) (22) 

with the replacement p — > p(t), x x{t), and {p(r), x(r)} is the particle 
trajectory in phase space. We will discuss the hermiticity property of the 
fractional Hamiltonian in Sec. IV. 

The exponential in Eq. ipTi) can be written as exp{^Sa{p,x)} if we intro- 
duce the fractional canonical action for the trajectory {p{t), x{t)} in phase 
space 

Sa{p,x) = J dT{p{r)xiT) - if„(p(r),x(r)). (23) 

ta 

Since the coordinates xq and xat in definition (!20l) are fixed at their initial 
and final points Xq = Xa and xn = Xb, all possible trajectories in Eq. (l23l) 
satisfy the boundary condition x(tb) = Xb and x(ta) = Xa- We see that the 
definition given by Eg. (120 p includes one more pj integral than Xj integral. 
Indeed, while xq and xn are held fixed and the Xj integrals are done for 
j = l,...,A^ — 1, each increment Xj — is accompanied by one pj integral 
for j = 1,...,N. The above observed asymmetry is a consequence of the 
particular boundary condition. That is, the end points are fixed in position 
(coordinate) space. There exists the possibility of proceeding in a conjugate 
way, keeping the initial pa and final pb momenta and fixed. The associated 
kernel can be derived going through the same steps as before, but working 
in the momentum representation (see, for example, Ref. [18]). 

Taking into account Eq. fll7p it is easy to check directly the consistency 
condition 

oo 

Kfl\xbtb\Xata) = j dx'Kf\xbtb\x't') ■ Kf\x't'\Xata)- 
— oo 

This is a special case of the general fractional quantum-mechanical rule: 
amplitudes for events occurring in succession in time multiply 

oo 

KL{Xbtb\Xata) = j dx' K L{Xbtb\x't') ■ K L{x't' \Xata) ■ (24) 
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4 Fractional Schrodinger equation 

The kernel KL{xbtb\xata) which is defined by Eqs.( !T2|) and (fT3!) . describes the 
evolution of the fractional quantum-mechanical system 

oo 

'l(jj-{Xb,tb) = j dXaKL{Xbtb\Xata) ■ 1pi{Xa,ta), (25) 

— oo 

where iljj^{xa,ta) is the fractional wave function of the initial (at the t = to) 
state, and ipf{xb, h) is the fractional wave function of the final (at the t = tb) 
state. 

In order to obtain the differential equation for the fractional wave function 
ipi^Xyt), we apply Eq. (!25!) in the special case that the time tb differs only by 
an infinitesimal interval e from ta 



'ijj{x,t + e)= / dyKLix,t + e\y,t) ■ i){y,t). 



t+T 

Using the Feynman's approximation J drV{x{T)) ^ £:F[(a; + y)/2] and 



t 



the definition given by Eq. flTTl) we have 



oo oo 

^{x,t+e)= j dy^ I dpeM^^-^^'^^^-yV{^y)}.^{y,t). 



We may expand the left- and the right-hand sides in power series 



oo oo 

^l,ix,t)+e^^= ! dy\ [ dpe^'-^ {1 - x (26) 



dt J 271% J ' n 

— oo — oo 

(i-^^n^))-^(y,t)- 

Then, taking into account the definitions of the Fourier transforms. 
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-oo 



oo oo 
— oo — oo 

and introducing the quantum Riesz fractional derivative (^V)" 

oo 

(W)Xx,t) = -^ j dpe^'^lplMP^t), (27) 
we obtain from Eq. fl26l) . 



^(x, t) + = ^(a;^ t) + i^{hVYiP{x, t) - ^-eV{x)^lj{x, t). 

This will be true to order e if ilj{x,t) satisfies the fractional differential 
equation 

ih^ = -D^inwy^ + V{x)^. (28) 
ot 

This is the fractional Schrodinger equation for a fractional quantum par- 
ticle moving in one dimension. 

Equation (l28l) may be rewritten in the operator form, namely 

ih^ = H^^, (29) 
where Ha is the fractional Hamiltonian operator: 

Ha = -DaihV)'' + V{x). (30) 

Using definition fl27|) one may rewrite the fractional Hamiltonian in 
the form given by Eq. fl22|) . 

The Hamiltonian Ha is the Hermitian operator in the space with scalar 
product 



X) = dx(f)*{x,t)x{x,t). 
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To prove the hermiticity of Ha, let us note that in accordance with the 
definition of the quantum Riesz fractional derivative given by Eg. fl27p there 
exists the integration-by parts formula 

((/.,(wrx) = ((w)>,x). (31) 

The average energy of fractional quantum system, with Hamiltonian H^, 

is 

oo 

Ec,= j dxip*{x,t)Haipix,t). (32) 

— oo 

Taking into account Eq. flMl) we have 



Ea= / dxij*{x,t)Haij{x,t) = / dx{H^2P{x,t)yij{x,t) = E*, 



and, as a physical consequence, the energy of a system is real. Thus the 
fractional Hamiltonian defined by Eq. (l30l) is the Hermitian or self-adjoint 
operator 

{H^<P,x) = {<P,HaX)- 

Since the kernel Ki^Xbtblxata), thought of as a function of variables Xb 
and tb, is a special wave function (for a particle which starts at Xa, ta), we see 
that Kl must also satisfy a fractional Schrodinger equation. Thus, for the 
quantum system described by the fractional Hamiltonian [Eq. (l30|l ]. we have 



d 

lh—KL{Xbtb\Xata) = - D a{flV bT K L{Xbtb\Xata) + V {Xb)K L{Xbtb\Xata) , h > ta, 



dt 



b 



where the low index "6" means that the quantum fractional derivative 
acts on the variable Xb- 
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5 Free particle. Fractional uncertainty rela- 
tion 



As a first physical application of the developed FQM and the fractional 
Schrodinger equation (1251) . let us consider a free particle. The fractional 
Schrodinger equation for a free particle has the fractional plane wave solu- 
tion 

ilj{x,t) = C -expU— -t '> , (33) 

where C is a normalization constant. In special Gaussian case (a = 2 and 
D2 = l/2m) Eq. fl33l) gives a plane wave of the standard quantum mechanics. 
Localized states are obtained by a superposition of plane waves 

00 

— 00 

Here (p{p) is the "weight" function. We will study Eq. (lMI) for a one- 
dimensional fractional Levy wave packet, 



00 



^d^^t) = 7^ / dpexp \ ^ . exp <^ z— - } , (35) 



2nh J ' { 2%" I \ n 

— CO 

with the "weight" function 

ip{p) = A„ ■ exp <^ — } , Po>0, ly <a, 



2r 

where A^, is a constant, / is a space scale and a is the Levy index, 1 < 
a<2. 

In the following we will be interested in the probability density p{x, t) 
that a particle occupies a position x, and the probability density w{p, t) that 
a particle has particular values p of the momentum. The wave function 
%lji^{x,t) defined by Eq. fl35l) . gives the probability density p{x,t) 

00 

p{x,t) = \^p,{x,t)\' = J^^ I rfMp2exp|-^l^|^|x (36) 
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Now, we can fix the factor A^, sucli tliat / dxp{x, t) = J dxlipi^^x, t)p = 1, 
witfi tlie result 



where r(^) is the gamma functioiJl]. The relationship between the prob- 
ability densities p{x,t) and w{p,t) may be derived from the relationship be- 
tween fractional wave functions in the space il)]^{x,t) and momentum (pipjt) 
representations, 

oo 
— oo 

where we have 

(f){p, t) = exp I — I ■ exp I -t — I . (39) 

Note that (j){p, t) satisfies the fractional free particle Schrodinger equation 
in the momentum representation 



ih — -- — = D^\p\ ^{p,t), 4>ip,0) = exp 



One then obtains 



oo oo 



dx\ipj^{x,t)\''^ = I dx J dpdp' exp i^i— — ^ '^'^ j- (f){p, t)(j)*{p', t) = 

(40) 



-oo — oo 



^The gamma function T{z) has the famihar integral representation r(z) — J dtV- "^t 



Rez > 0. 
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{27rh) 

because of 



dp\c^{p,t)\' 



{2nh) 



oo 

J dxexpi^i— — j> = 6{p — p'). 



Equation ([30]) suggests, for the probability density in momentum space, 
the following definition: 

w{p,t) = f-\<f>{p,t)\'. (41) 
Thus, for the momentum probability density w{p,t), we have 

w{p, t) = w{p) = -TT^TYT exp <^ — } . (42) 



2^r(i) - I h ^ 

This is time independent, since we are considering a free particle. 

In coordinate space the probability of finding a particle at the position 
X in the "box" dx is given by p{x,t)dx. Correspondingly, the probability of 
finding the particle with momentum p in dp is represented by w{p,t)dp. 

We are also interested in the average values and the mean-// deviations 
of position and momentum for the present probability densities defined by 
Eqs.( l36|) and (l42l) . The expectation value of the space position can be calcu- 
lated as 



< X >= / dxxp{x,t) = (43) 



^ dxx I dpdp e:!Cp\'^~ — 0(P) ^)0*(p') ^)- 



— oo — oo 

Making the substitution 

h d 
I op 
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we will have 



oo oo 

<^>-(2ip / /**'(f|exp{,(^}),(p,«),-(p'.0^ 

— OO — oo 

Integrating by parts gives 



oo 



<x >= - - / dp[— — \p-po\ ttIpI exp ' 



— oo 

It is easy to check that the first term in the brackets vanishes, and we 
find that the position expectation value is 

<x>= aD^po'^-H. (44) 
Using the dispersion relation given by Eg. (1191) . we may rewrite < x > as 

<x >= ■ t = vot. (45) 

Here vq = {dEp/dp)\p=pg is the group velocity of the wave packet. We see 
that the maximum of the Levy wave packet [(Eq j35l) ] moves with the group 
velocity vq like a classical particle. 

The mean-/i deviations (/i < z/) of space position < |Aa;|^ > is defined by 

oo 

<iA.r>^<i.-<.>r>^/..i.-<.>rp(.,t)^ 



- — — -r / dx\x— < X > r / apap exp < i 

{^-nny J J V n 



This equation can be rewritten as 



0(p,t)0*(p',t). 



< \ Ax\^ >= -Af{a, /i, u; T, r^o), (46) 
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where we introduce the following notations 



oo oo 



l^l^y f f f 

Af{a,fi,u]T,7]Q) = ^-^^ J rf^l^l'' y dr] J rfr/'exp{i(r7-V)(?+arr/o°-^)}x 

— oo — oo — oo 

(47) 

exp{-zr(|r^r " WH - 1^ - ^ol' " W - ^oH 



Pol f2^^^h 



and 



2y-h' ^ h \ I 

So, for the /z-root of the mean-// deviation of position (space position 
uncertainty for the Levy wave packet), we find 

< \Axr >'/^'= ^Af'/^{a,f,,u;T,^,). (48) 

Further, with Eq.( l42l) the expectation value of the momentum is calcu- 
lated as 



< P >= / dppw{p) = / dp{p — po)w{p) + / dppow{p). (49) 



The first integral vanishes, since w{p) is an even function of (p — po), and 
the momentum expectation value is 

<P>= Po- (50) 
The mean-/x deviation of the momentum is 

< \Apr >= j dp\p-<p> rw{p) = l^jj -j^. (51) 

— oo ^ 

Then the momentum uncertainty (the /x-root of the mean-// deviation of 
momentum) is 
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Together with Eq.fHHl). this leads to 




(52) 



< lAxI'^ >^'^'< lApl'^ > 



2V^ 



U^/^^{a,iJi,u-T,r^^), (53) 



< u < a. 

where A/'(a, /x, z/; r, ?7q) is given by Eq. (l47|l . 

This relation implies that a spatially extended Levy (or fractional) wave 
packet corresponds to a narrow momentum spectrum, whereas sharp Levy 
wave packet corresponds to a broad momentum spectrum. 

Since Af{a,n,u; t,i]q) > 1 and r(^)/r(i) ^ 1/z/ Eq.([53]) becomes 



< \Axr |Aj9r /i < «, 1< « < 2, (54) 

with V = a. 

Note that for the special case when a = 2 we can set = a = 2. 
Thus, for the standard quantum mechanics (a = 2) with the definition of 
the uncertainty as the square-root of the mean-square deviation, Eq. flS^ was 
established by Heisenberg [19], (see, for instance, ref.[20]). The uncertainty 
relation given by Eq.f l54l) can be considered as fractional generalization of the 
well known Heisenberg uncertainty relation. Thus Eqs.f ll2l) -( |T5|) . (l2Tl) - (!2^ . 
( !28|) . (!30l) and (!54l) are the basic equations for the new FQM. 



6 Fractional statistical mechanics 

In order to develop the fractional statistical mechanics (FSM), let us go 
in the previous quantum-mechanical consideration from imaginary time to 
"inverse temperature" (3 = l/fc^T, where /c^ is Boltzmann's constant and T 
is the temperature, it —>■ hj3. In the framework of the traditional functional 
approach to the statistical mechanics, we have the functional over the Wiener 
measure [5], [IB], [21], 
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x(/3)=x hp 
p{x,P\Xo)= j VwienerXiu) ■ exp{-^ j duV{x{u))}, (55) 
x{0)=xo 

where p{x,P\xo) is the density matrix of the statistical system in the 
external field V{x), and the Wiener measure [6] generated by the process of 
the Brownian motion is given by 

\ V WienerX [U )...= \im^ X (56) 



x{0)=xo 

N 

r I ^ 

[Xj — Xj-i 



X JJ^exp 



m , .2 



here ^ = h[3/N and x^q = x. 

The FSM deals with Levy or fractional density matrix p^(x, /?|xo) which 
is defined by 



r:(pj=x np 

j Pl^^j^x(u) ■ exp{-i j duV{x{u))}, (57) 



x{i3)=x np 

Pl{x,(3\xq) = 

x(0)=xo 

where we introduce the new fractional functional measure (we will call 
this measure by the Levy functional measure) by 

x{P)=x 

j VLevyx{u). ..= \im j dxi...dxN-i{fi''~^D^q)-^''^x (58) 



x(0)=xo 

N 



" 1 



here ^ = hP/N, xn = x and the Levy function is given by Eq.([T] 
Equations fl57j) and fl58l) define the fractional quantum statistics via new Levy 
path integral. 

The partition function Z or free energy F, Z = e~^^ is expressed as a 
trace of the density matrix pi^{x, P\xo): 
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Z = e-^^ = j dxpL{x,p\x) = 
j dx j Visvyxiu) ■ exp{-^ j duV{x{u))}. 

x{0)=x{(3)=x 

With the definition fl2U]) the equation for the partition function becomes 

Z = e-^^ = (59) 



hl3 

J J Dx{t) j Dp(r) exp{ — — J du ^^—ip{u)x{u) + Ha{p{u), x{u))^ , 

x{0)=x{l3)=x 

where the fractional Hamiltonian H^ip, x) has form of Eq. fl2^ . and p{u),x{u) 
may be considered as paths running along on "imaginary time axis", u = it. 
The exponential expression of Eq. (!59|) is very similar to the fractional canon- 
ical action [Eq. (l23l) ]. Since it governs the fractional quantum-statistical path 
integrals, it may be called the fractional quantum-statistical action or frac- 
tional Euclidean action, indicated (following Ref . [T5] ) by the superscript (e). 



The parameter u is not the true time in any sense. It is just a parameter 
in an expression for the density matrix (see, for instance, Ref. [3]). Let us 
call u the "time", leaving the quotation marks to remind us that it is not 
real time (although u does have the dimension of time). Likewise x{u) will 
be called the "coordinate" and p{u) the "momentum". Then Eq. fl57l) may 
be interpreted in following way. 

Consider all possible paths by which the system can travel between the 
initial x(0) and final x{(3) configurations in the "time" hjS. The fractional 
density matrix is a path integral over all possible paths, the contribution 
from a particular path being the "time" integral of the canonical action 
[considered as the functional of the path {p{u),x{u)} in the phase space] 
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divided by fi. The partition function is derived by integrating over only 
those paths for which initial x(0) and final x{f3) configurations are the same, 
and after that we integrate over all possible initial (or final) configurations. 

The fractional density matrix p^^(x, /5|xo) of a free particle {V = 0) can 
be written as 

oo 

pP{x,(3\xo) = ^ I dpexp|z^^^^-/3D„br|= (60) 

) 

\x-xo\ (1, 1/a), (1,1/2)" 

HD^py/- ' (1,1), (1,1/2) J ' 

where is the Fox's H function (see, Refs. [T5] - \T7]). 
For a linear system of space scale Q the trace of Eq. (!60l) leads to 



a\x — xo\ 



1,1 
2,2 



00 



n —00 

When a = 2 and D2 = l/2m. Eg. (1601) gives the well-known density matrix 
for a one- dimensional free particle (see Eq. (10-46) of Ref. [3j or Eq.(2-61) of 
Ref. [211): 



The Fourier representation P^l\p, l3\p') of the fractional density matrix 
pf\x, (3\xo) defined by 



pL\p^f^\p')= J dxdxop''°\x,(3\xo)exp{-^{px -p'xo)} 
—00 

can be rewritten as 

Pl\p,(^\p') = 2nh6{p~p') ■ e^'^^-IPl". 

In order to obtain a formula for the fractional partition function in the 
limit of fractional classical mechanics, let us study the case when hjS is small. 
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Repeating consider, similar to Feynman's (see, Chap. 10 of Ref. [5]) for the 
fractional density matrix pi^{x, f3\xQ) we can write the equation 



— oo 

Then the partition function in the limit of classical mechanics becomes 

oo ^ oo 

Z = / = 1 'i--^""'- (62) 

— OO — oo 

This simple form for the partition function is only an approximation, 
valid if the particles of the system cannot wander very far from their initial 
positions in the "time" h(3. The limit on the distance which the particles can 
wander before the approximation breaks down can be estimated in Eq. (l60ll . 
We see that if the final point differs from the initial point by as mush as 

Ax ~ hipD^y^'^ = n 

the exponential function of Eq. fl60l) becomes greatly reduced. From this, 
we can infer that intermediate points on paths which do not contribute 
greatly to the path integral of Eq. fl60l) . If the potential V{x) does not al- 
ter very much as x moves over this distance, then the fractional classical 
statistical mechanics is valid. 

The density matrix p^{x, j3\xo) obeys the fractional differential equation 




dp^ix, (3\xo)_ ^ _j^^^^^^yp^^^^ ^1^^) + V{x)pj^{x, P\xo) (63) 



d[3 
or 

9pj^(x,/?|xo) 



HapL{^,(^\Xo), Pi(x, 0|Xo) = 5{x-Xq), 



dp 

where the fractional Hamiltonian Ha is defined by Eq. fl30l) . 

Thus Eqs. ( l57l) - ( !60l) and ( !63l) are the basic equations for our FSM. 
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7 Conclusion 



We have developed a path integral approach to FQM and FSM. The approach 
is based on functional measures generated by the stochastic process of the 
Levy flights whose path fractional dimension is different from the fractional 
dimension of the Brownian paths. As was shown by Feynman and Hibbs, 
fractality (difference between topological and fractional dimensions) of the 
Brownian paths lead to standard (nonfractional) quantum mechanics and 
statistics. The fractality of the Levy paths as shown in the present paper 
leads to fractional quantum mechanics and statistics. Thus we develop a 
fractional background which leads to fractional (nonstandard) quantum and 
statistical mechanics. 

The Feynman quantum-mechanical and Wiener statistical mechanical 
path integrals are generalized, and as a result we have fractional quantum- 
mechanical and fractional statistical mechanical path integrals, respectively. 
A fractional generalization of the Schrodinger equation has been derived 
using the definition of the quantum Riesz fractional derivatives. We have 
defined the fractional Hamilton operator and proved its hermiticity. The 
relation between the energy and the momentum of nonrelativistic fractional 
quantum-mechanical particle has been found. The equation for the fractional 
plane wave function was obtained. We have derived a free particle quantum- 
mechanical kernel using Fox's H function. In the particular Gaussian case 
(a = 2), the fractional kernel takes the form of Feynman's well-known ker- 
nel. For the Levy wave packet the position and momentum uncertainties 
were calculated analytically. The fractional generalization of the Heisenberg 
uncertainty relation has been established. 

Equations. ([12)- (USD, (12I])-(I21, ([2H]), ([30]) and ([MD are the basic equations 
for our FQM. 

Following the general rule and replacing it hj3, we obtain the path 
integral formulation of the FSM. An equation for the fractional partition 
function has been derived, and the fractional quantum-statistical action in- 
troduced into the quantum statistical mechanics. The density matrix of a 
free particle has been expressed analytically in terms of the Fox's H function. 
It is shown that Eq. (l60|) for the fractional density matrix in a special Gaus- 
sian case (a = 2) gives the well-known equation for the density matrix of 
free particle in one dimension (see Eq.(2-61) of Ref. |21j). We have found the 
formula for the fractional partition function in the limit of fractional classical 
mechanics and discuss the validity of this formula. A fractional differential 
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equation of motion of density matrix has been established. Equations. (1571) 
- (1^ and are the basic equations for our FSM. 

We finally mention that the developed approach to quantum and statis- 
tical mechanics can easily be generalized to a ci- dimensional consideration, 
using, a d-dimensional generalization of the fractional and the Levy path 
integral measures. 
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